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Resume 



Lors de cette these mes recherches se sont principalement orientees vers 
I'etude de systemes antiferromagnetiques en presence de configurations antifer- 
ronnagnetiques exotiques caracterisees par un invariant topologique. 

De nnultiples tentatives infructueuses pour mener a terme le passage a une 
theorie des channps continue du modele antiferromagnetique de Heisenberg sur 
le reseau triangulaire dans la perspective de mettre en evidence la presence du 
ternne de Hopf dans le Lagrangien m'ont conduits a esquisser pour les reseaux 
antiferronnagnetiques de basse dimension la construction de nouveaux degres 
de liberte antiferromagnetiques locaux. Cette construction originale reproduit la 
hierarchie des algebres reelles a divisions verifiee par le modele o non-lineaire 
(limite continue du Hamiltonien de Heisenberg classique en dimension deux) et 
revele la presence d'un champs de jauge locale inedit. Surtout cette construc- 
tion algebrique m'amene a suggerer un nouveau processus microscopique qui 
serait a I'origine du terme de Hopf, d'etats antiferromagnetiques topologiques 
(caracterises par I'indice de Pontrjagin) et de la transmutation de Fermi-Bose 
(supraconductivite). 

En outre j'ai mis en evidence pour des membranes magnetiques souples de 
genre topologique spherique et torique en presence d'un soliton magnetique I'ex- 
istence d'un retrecissement global renfle localement aux vagues du soliton. L'o- 
rigine geometrique de ce phenomene inedit m'a amene a interpreter la frustration 
geometrique des membranes magnetiques souples comme la competition entre 
les deux ordres topologiques presents (decrits par I'indice de Pontrjagin du soliton 
et le genre de la membrane). 

Le processus microscopique et la competition topologique suggeres elargissent 
le champs d'investigation initial en fournissant une origine microscopique a la 
transmutation de Fermi-Bose et un moyen de traquer cette physique exotique. 



Abstract 



In this thesis I present my research on the exotic configurations of antiferromag- 
netic systems characterised by a topological invariant. 

It is expected that taking the continuum limit of the antiferromagnetic Heisen- 
berg model on the triangular lattice gives rise to a Lagrangian with a Hopf term. 
The research presented in this thesis shows that the continuum limit is much 
more difFicult to achieve than may be expected. This difFiculty led me to outline 
the construction of novel local antiferromagnetic degrees of freedom for low di- 
mensional antiferromagnetic lattices. This new construction reproduces the real 
division algebra hierarchy satisfied by the nonlinear a-model (which corresponds 
to the continuum limit of the Heisenberg Hamiltonian for classical isotropic anti- 
ferromagnets in two dimensions) and reveals the presence of a novel local gauge 
field. This algebraic construction leads me to suggest a new microscopic mech- 
anism which should give rise to the Hopf term, as well as antiferromagnetic 
topological states (characterised by the Pontrjagin index), and the Fermi-Bose 
transmutation (superconductivity). 

In this thesis I have also studied elastic magnetic vesicles of spherical and 
toric genus in the presence of a magnetic soliton. My studies reveal a global 
shrinking, with local swellings in the regions where the soliton presents a spin- 
flip. The geometrical origin of this novel phenomena led me to interpret the 
geometric frustration of magnetic vesicles as the competition between the two 
topological orders present (one described by the Pontrjagin index of the soliton 
and the other by the genus of the vesicle). 

The microscopic mechanism and topological competition suggested above 
go beyond the scope of this thesis, providing a microscopic explanation for the 
Fermi-Bose transmutation and a way to deal with this type of exotic physics. 
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Introduction 



Le mecanisme responsable de la supraconductivite a haute temperature n'est 
pas encore connu : les theories « conventionnelles » de la matiere condensee 
ne predisent pas de temperatures critiques appartenant au domaine actuelle- 
ment atteint par les experiences (de AOK a 120K). A basse temperature I'in- 
teraction entre les electrons (les porteurs de charge du reseau cristallin) et les 
phonons (les vibrations elementaires du reseau cristallin) induit la supraconduc- 
tivite « ordinaire ». Par contraste I'etat supraconducteur a haute temperature 
serait du aux spins des electrons. 

Dans ce contexte un argument simple de topologie suggere I'existence pour un 
reseau bidimensionnel antiferromagnetique d'etats caracterises par un invariant 
topologique : il s'agit en fait de recouvrir la sphere unite (la variete des spins) 
en parcourant I'hypersphere unite (le plan-temps compactifie) — le nombre 
de recouvrement etant invariant topologique de Hopf. Une theorie de champs 
continue d'un tel systeme est generalement decrite par le modele a non-lineaire. 
I'emergence lors du passage a une theorie de champs continue de invariant de 
Hopf dans le Lagrangien (conjecturee par Dzyaloshinski, Polyakov et Wiegmann) 
est I'objet de multiples controverses : cela modifierait la statistique du systeme 
et revelerait I'existence d'une transmutation de Fermi-Bose [P0I88]. Pour les sys- 
temes de basse dimension ce mecanisme se generaliserait en suivant la hierarchie 
des algebres a divisions [Tze88]. Avant d'envisager tout passage a une theorie 
des champs continue, la construction d'un parametre d'ordre local pertinent est 
primordiale : la premiere partie de ce memoire suggere la construction d'un nou- 
veau parametre d'ordre local qui reproduit la hierarchie des algebres a divisions 
et dont le Lagrangien peut exhiber le terme de Hopf. 

De plus les etats statiques d'un reseau bidimensionnel antiferromagnetique 
sont ici caracterise par I'indice de Pontrjagin qui decompte les recouvrements de 
la variete des spins lorsque le reseau bidimensionnel est parcouru : a la limite con- 
tinue il s'agit d'envoyer une membrane compactifiee (la surface magnetique) sur 
la sphere unite S'^ (la variete des spins) via le modele a non-lineaire (le Hamilto- 
nien magnetique). Certaines membranes courbes admettent des distributions de 
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spins de type soliton caracterisees par un parametre de forme. La deuxieme partie 
de ce memoire etudie la competition entre I'energie magnetique de la distribu- 
tion des spins et I'energie elastique de la membrane en suggerant une approche 
topologique. Enfin un retrecissement global avec des renflements localises aux 
solitons est mis en evidence. 



Chapitre 1 

Le reseau antiferromagnetique 
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Le reseau antiferromagnetique 



1 Le modele antiferromagnetique de Heisenberg 

Dans les solides les interactions entre electrons sont souvent importantes et 
d'une extraordinaire complexite. En physique de la matiere condensee, ce prob- 
lems d'interaction electronique se reduit dans la plupart des cas a un problems 
ds spins coupiss sur rsssau [FraQl, Aus94, L6v97]. Ici les spins seront identiques 
et occuperont les noeuds de reseaux reguliers isotropes de basse dimension. 

1.1 Spins sur reseau 

L'etude du magnetisme necessite rintroduction d'un objet physique typ- 
iquement quantique : le moment angulaire intrinseque communement appele 
spin. Du point de vue de la theorie des champs, le spin d'une particule elemen- 
taire correspond a la « charge » de Noether associee a I'isotropie de I'espace- 
temps propre a cette particule elementaire [BL93, G0I8O, Siv95]. Autrement 
dit, le spin d'une particule elementaire traduit la rotation de ladite particule 
elementaire sur elle-meme, sa rotation dans I'espace-temps etant decrite par 
son moment angulaire orbital. Notons au passage que le moment angulaire 
orbital possede un equivalent classique contrairement au moment angulaire 
intrinseque. Dans le cadre present, a chaque noeud (i) du reseau est associe 
un spin represente par un operateur vectoriel Sj'* qui engendre une algebre de 
Lie [SW86] de dimension s selon la relation de commutation 

^\He^^^5ijSi^ (1.1) 

avec h la constante de Planck, c"'^" le tenseur unitaire parfaitement anti- 
symetrique d'ordre trois et 5ij le symbole de Kronecker. De fait V operateur 
de Casimir [SW86] fixe la dimension du degres de liberie local du reseau : 

S' = S,"Si, = s(s + l) Vi (1.2) 

Les reseaux envisages sont les pavages compacts regulier (sans lacune ni 
recouvrement) des espaces Euclidiens de basse dimension (d = 1, 2, 3, 4) avec 
des polytopes reguliers identiques [Cox73, Siv95]. Dans le cadre present, les 
reseaux seront done caracterises par un polytope regulier : 

- le reseau aura la dimension d du polytope regulier ; 

- les symetries qui preservent le polytope regulier preserveront le reseau ; 

- les translations qui reproduisent le reseau a partir d'un exemplaire du 
polytope regulier preserveront le reseau ; 

- la longueur des aretes du polytope regulier correspondra au pas a du 
reseau. 
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Observant que les polytopes reguliers de dimension d sont egalement des 
pavages reguliers de la sphere S^~^ [dV64, Cox73], introduisons la notion de 
« sur-reseau » : nous entendrons par « sur-reseau » de dimension d un reseau 
de meme dimension dont pent s'extraire un reseau regulier de dimension d en 
assimilant ses « sur-nceuds » a des spheres S^~^ identiquement et reguliere- 
ment pavees (Figure 1.1). II est evident qu'a un reseau regulier pent etre 
associe plusieurs « sur-reseaux » , et reciproquement. 

Le comportement magnetique sur de longues distances (longues devant 
le pas a du reseau) sera apprehende en passant a une theorie de champs 
continue selon une prescription due a Affleck [Aff88, Alf89]. II s'agit de faire 
tendre conjointement le pas a du reseau vers zero et le nombre quantique s 
vers I'infini tout en maintenant constantes les quantites physiques mesurables 
experimentalement. Cette prescription traite implicitement les spins a la lim- 
ite dassique. En efFet la relation de commutation (1.1), qui doit s'ecrire pour 
garder un sens 



0' 



'J 



i-e"^^5,,-^, (1.3) 
s s 



s'evanouit lorsque le nombre quantique s s'approche de I'infini. Nous sup- 
poserons les resultats physiques obtenus applicables pour des valeurs physiques 
du nombre quantique s (s = ^ ou s = 1) et des valeurs du pas a carac- 
teristiques des solides cristallins (de I'ordre de 4 A). II est evident que la 
construction d'un parametre d'ordre local pertinent est un prelude necessaire 
et complexe qui doit a la fois clarifier les phenomenes physiques observes et 
valider la justesse du passage a une theorie de champs continue. 



1.2 Interactions d'echange : Hamiltonien de Heisenberg 

Le magnetisme electronique provient de V interaction Coulombienne en- 
tre electrons qui force les spins dans des etats ordonnes : le principe de 
Pauli impose des etats completement antisymetriques lorsque deux femiions 
sont echangees. Generalement le Hamiltonien initiale d'un solide cristallin 
qui est extraordinairement complexe pent etre reduit sous une forme simple 
parametrisee seulement en termes de spins. Le plus souvent cet Hamiltonien 
effectif est de Heisenberg [Lev97] : 

Mh = - 51 Jif SiaSj^, (1-4) 

ou les indices i et j etiquettent les porteurs de spins. Les tenseurs generiques 
J^^ decrivent les interactions d'echange. Pour un couplage antiferromag- 
netique isotrope entre premiers voisins sur un reseau regulier de spins nous 
lisons 

=^if=/ES^°S,-a, (1.5) 

{h3) 

avec J une constante de couplage spin-spin positive. 
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A haute temperature (temperature superieure a la temperature de Curie) 
il n'y a pas d' ordonnancement magnetique : le systeme magnetique est dans 
un etat paramagnetique. Le systeme de spins est alors localement isotrope et 
invariant par rotation : un systeme paramagnetique possede les symetries de 
son reseau de spins. A basse temperature differents ordres peuvent emerger. 
A tres basse temperature (temperature proche du zero ahsolu) les systemes 
de spins reguliers et isotropes acquierent une orientation particuliere dans 
I'espace : il y a un ordonnancement magnetique. Cette direction spontanee 
brise I'invariance locale par rotation du systeme. En observant qu'un sys- 
teme physique ne pent pas changer graduellement de symetrie (une symetrie 
existe on n'existe pas) Landau a mis en evidence qu'une transition de phase 
du second ordre doit separer les etats de differentes symetries [Lev97, Siv95]. 
Generalement quand une symetrie est brisee Vetat non-symetrique est carac- 
terise par un parametre d 'ordre local : selon Landau tout parametre local dont 
la valeur moyenne est nulle dans I'etat symetrique et non-nulle dans I'etat 
non-symetrique pent constituer un parametre d 'ordre local [Siv95]. 




Le reseau triangulaire paramagnetique Le reseau carre paramagnetique 



Figure 1.1: Les reseaux triangulaire et carre paramagnetiques : I'invari- 
ance locale par rotation des deux systemes magnetiques est schematisee en 
magnifiant les noeuds porteurs de spin. Pour chaque reseau le « sur-pavage » 
represente un choix possible de « sur-reseau » : les plaquettes grisees corre- 
spondant alors aux « sur-noeuds » . 
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1.3 Etat antiferromagnetique : sous-reseaux de Neel 

Dans un solide cristallin antiferromagnetique deux spins adjacents font 
de leur mieux pour s'aligner selon des orientations opposees (||). II se cree 
ainsi un motif alterne a travers le reseau forme par les spins. Sur un reseau 
regulier de spins cet arrangement alterne se traduit par I'emergence d'une an- 
tisymetrie locale [Siv95] : les noeuds de chaque polytope d'un «sur-reseau» 
du reseau regulier initial sont colores. Cette brisure de symetrie locale carac- 
terise Vetat antiferromagnetique. Pour illustration : la chame paramagnetique 
de pas a (Figure 1.2) est preservee par la translation de pas a, alors que 
la chame antiferromagnetique de pas a (Figure 1.3) est alternee par I'anti- 
translation de pas a et preservee par la translation de pas 2a (i.e. le pas des 
« sur-chaines »). L'etat classique correspondant a I'energie la plus basse du 
Hamiltonien de Heisenberg antiferromagnetique (1.5), communement appele 
Vetat de Neel, constitue une representation immediate du reseau de spins col- 
ore. Les spins du reseau de Neel ayant la meme orientation forment les sous- 
reseaux de Neel : les reseaux cubiques^ de Neel sont bipartites (Figures 1.3 
et 1.5), le reseau triangulaire de Neel est tripartite (Figure 1.4). Enfin les 
«sur-noeuds» du reseau de Neel sont usuellement appeles les plaquettes de 
Neel : chaque plaquette de Neel correspond en general a une description lo- 
cale des sous-reseaux de Neel, d'ou une certaine ambiguite qui sera levee en 
utilisant systematiquement la notion de « sur-noeud » dont le caractere local 
est sans equivoque (Figures 1.1 et 1.2). 



• — I — • • — I — • • — I — • • — I — • • — I — • • — I — • • — I — • • — I — • 



Figure 1.2: La chaine paramagnetique. 



t ' i t ' i M - I'i f4H-4-|-'-i- - 1-4 

La chaine de Neel 



i ' l i ' t i ' t ^-M-'^^4-4^4-4-'4 

L'anti-chame de Neel 
Figure 1.3: La chaine ferromagnetique. 



^Dans le cadre present : la chame (d=l), le reseau carre (d=2), le reseau cubique (d=:3) 
et le reseau 4-cubique (d=4). 
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Le reseau triangulaire de Neel 



L'anti-reseau triangulaire de Neel 



Figure 1.4: Le reseau carre ferromagnetique : rarrangement alterne des 
spins est tricolore. L'antisymetrie tricoloree est usuellement representee par 
trois families de vecteurs dont les directions prises deux a deux forment un 
angle de ^vr : les sous-reseaux de Neel (— \ et /). 



I ($^^ I ($^^ I ($^^ I ^ I I I iSl^ ^-h^ 





















Le reseau carre de Neel 



L'anti-reseau carre de Neel 



Figure 1.5: Le reseau carre ferromagnetique : Tarrangement alterne des 
spins est bicolore. L'antisymetrie bicoloree est usuellement representee par 
deux families de vecteurs de direction identique mais de sens oppose : les 
sous-reseaux de Neel (© et (^). 



1 Le modele antiferromagnetique de Heisenberg 



9 



Par construction deux antioperations successives laissent le reseau de 
spins invariant : nous caracteriserons les « sur-reseaux » par leur helicite 
[Siv95]. Consequences immediates : toute antioperation autre que I'anti- 
translation n'affecte que les «sur-noeuds» tout en changeant V helicite du 
« sur-reseau », et deux antitranslations successives ne constituent rien d'autre 
qu'une translation du « sur-reseau ». En d'autres termes, puisque I'antisymetrie 
caracterise I'etat antiferromagnetique, les « sur-reseaux » sont transparents a 
Fantiferromagnetisme du systeme (a I'lielicite du « sur-reseau » pres) alors 
que Fantiferromagnetisme est concentre en leurs « sur-noeuds » . A posteriori 
le « sur-reseau » apparait done comme le reseau de I'etat antiferromagnetique : 
un parametre d'ordre local pertinent de Vetat antiferromagnetique doit etre 
associe a un «sur-nceud» [DR88, DR89, ADDJ93]. 

1.4 Antiferromagnetisme frustre et assouvi 

Deux grands types de frustrations existent pour les systemes antiferro- 
magnetiques [Kaw98, Lev97] : 

- il existe une competition entre plusieurs ordres magnetiques de dif- 
ferentes natures dont un au moins est antiferromagnetique : e.g. un 
ordre antiferromagnetique entre premiers voisins s'opposant a une in- 
teraction d'echange ferromagnetique entre seconds voisins; 

- le reseau interdit la satisfaction de toutes les interactions d'echange : 
e.g. sur le reseau triangulaire il est impossible d'aligner antiparallele- 
ment trois spins deux a deux adjacents. 

Seule la frustration due au reseau est pertinente dans le contexte present. 
Cette frustration geometrique est caracterise par im nombre de sous-reseaux 
de Neel superieur a deux. Le « sur-reseau » est done transparent a la frustra- 
tion geometrique [DR89, ADDJ93]. 
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2 Le « sur-reseau » antiferromagnetique 

Dans la partie precedente le reseau des porteurs de spin apparaTt comme 
le reseau de I'etat paramagnetique, le « sur-reseau » comme le reseau de I'etat 
antiferromagnetique — les «sur-nceuds» renfermant la non-symetrie. 

J'esquisse une mise en evidence pour les « sur-reseaux » antiferromagnetiques 
reguliers de basse dimension d une hierarchie algebrique analogue a celle du mod- 
ele a non-lineaire [Tze88, TN89, GT96] : les «sur-noeuds» se comportent ici 
comme des « sur-particules » a part-entiere, les « sur-reseaux » antiferromagne- 
tiques comme des reseaux de « sur-particules ». 

2.1 Les degres de liberte antiferromagnetiques locaux 

Je construis les degres de liberte antiferromagnetiques locaux associes a 
un «sur-noeud» comme etant les combinaisons lineaires des degres de lib- 
erie locaux {i.e. les spins) associes au polytope regulier du «sur-noeud» 
qui reduisent les representations irreductibles du groupe des transforma- 
tions laissant le polytope regulier du «sur-noeud» invariant. II s'avere alors 
avantageux de se representer le polytope regulier du «sur-noeud» comme 
un pavage regulier de la sphere S'^^^ avec d la dimension du reseau. Or 
chaque groupe des transformations laissant invariant un pavage regulier de 
la sphere S'^^^ admet une representation dans Valgebre reelle a divisions^ 
de dimension 2^'^ [EHH+90, GT96, Dix94, dV64, Cox91] : les combinaisons 
lineaire seront construites dans cette algebre A^. Je choisi a priori comme 
parametre d'ordre antiferromagnetique local la combinaison lineaire algebrique 
normee a la limite continue qui est preservee par I'antioperation caracterisant 
I'etat antiferromagnetique a I'helicite pres. Lors du passage a une theorie 
de champs continue je suppose que les combinaisons lineaires algebriques 
restantes se comportent comme si les degres de liberte {i.e. les spins) du 
polytope regulier etaient effectivement decrits par un champs continu. Enfin 
j 'impose une representation locale : les degres de liberte antiferromagnetiques 
locaux seront definis a une rotation locale pres. La rotation locale s'appliquant 
aux degres de liberte du «sur-nceud» et non a sa sphere S^"^ doit s'inter- 
preter comme un champ de jauge locale [Nab97, BL93, Nak90]. Cette derniere 
hypothese est la plus originale : la representation implicitement choisie est 
generalement globale [DR88, DR89, ADD.J93]. 

Les interactions d'echange n'engendrent pas ici un ordonnancement an- 
tiferromagnetique global (c'est-a-dire «a la Neel») mais forcent les spins a 
s'organiser en « sur-nceuds » : les spins appartenant a un meme « sur-nceud » 

^Real Division Algebra 



2 Le « sur-reseau » antiferromagnetique 
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se lient entre eux tout en se desolidarisant des spins appartenant aux « sur- 
noeuds » voisins. Dans ce contexte les « sur-noeuds » acquierent une legitimite 
physique immediate : ils se comportent comme des « sur-particules ». Dans ce 
nouvel etat antiferromagnetique le comportement magnetique sur de longues 
distances (longues devant le rayon d'une « sur-particule ») provient des in- 
teractions entre les « sur-particules » decrites ici comme une desolidarisation 
entre spins lies a deux « sur-particules » adjacentes. En outre, un polytope 
regulier de dimension d n'etant rien d'autre qu'un reseau regulier sur la 
sphere S^~^, une « sur-particule » n'est rien d'autre qu'un reseau antiferro- 
magnetique regulier de la sphere S^~^ : I'utilisation de la theorie des groupes 
pour effectuer un changement de degres de liberte apparait done legitime. 
Une representation etant definie a une rotation pres, I'antisymetrie antiferro- 
magnetique de chaque reseau spherique {i.e. de chaque « sur-particule ») est 
definie a une rotation pres : les reseaux spheriques sont desolidarises en im- 
posant un champs de jauge locale. Les rotations invoquees sont celles qui lais- 
sent la sphere S''^"^ invariante. La theorie des groupes permet ainsi de separer 
Vantisymetrie imposee par les interactions d'echange et la symetrie heritee du 
polytope regulier, le reseau maintenant une certaine cohesion traduite par 
I'existence d'un champs de jauge locale. Aussi les degres de liberte antifer- 
romagnetiques locaux heritent-ils de la hierarchic algebrique des groupes des 
transformations preservant les polytopes reguliers, c'est-a-dire la hierarchic 
algebrique des algebras reelles a divisions A^. Nous nous attendons done a 
rencontrer une physique exotique. 

2.2 Cas unidimensionnel : la chaine de Haldane 

La chaine reguliere est pavee avec le 1-cube C° represente dans V algehre 
reelle a divisions M par le doublet 



{+1,-1}, 



(1.6a) 
(1.6b) 




Figure 1.6: Representation reelle du 1-cube C°. 
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Les transformations preservant le 1-cube se reduisent a I'identite 1 et a 
la rotation/antisymetrie —1. Le choix d'une «sur-chaine» etant fait, pour 
chaque «sur-noeud» (p) de la « sur-chaine » nous etiquetterons par (po) et 
(pi) les deux noeuds representes respectivement par les reels +1 et —1. Re- 
spectant la prescription introduite dans la Sous-Section 2.1, nous introduisons 
le parametre d'ordre antiferromagnetique local rip" tel que 

2i/n/ = Sp„"-Sp,", (1.7a) 
et la representation triviale Ip" telle que 

2 01^^^ = V + S^,^ (1.7b) 

A la limite classique I'operateur de Casimir (1.2) impose entre les nou- 
veaux degres de liberte antiferromagnetiques locaux les relations exactes suiv- 
antes : 

nl=l-^ll, (1.8a) 

= 0' (l-8b) 

avec 

u^^s{s + l). (1.9) 

Les deux contraintes (1.8) montrent clairement que le nombre total de de- 
gres de liberte (egal a quatre) est conserve lors du changement de degres de 
liberte (1.7). En outre lorsque le pas a du reseau et le nombre quantique s 
tendent conjointement vers zero et I'infini respectivement (prescription d'Af- 
fleck [Aff88, Aff89]) le parametre d'ordre antiferromagnetique local tip" est 
effectivement norme. Observons enfin que pour la theorie des champs conti- 
nus ainsi construite le parametre d'ordre local vit sur la sphere unite ^S"^. 

Etant donnes les degres de liberte antiferromagnetiques locaux Tip" et 
Ip" decrivant ime « snr-particule » (p) d'une chaine antiferromagnetique, les 
spins internes de ladite « sur-particule » Spg" et Sp^" s'obtiennent aisement en 
inversant (1.7) ; nous lisons 

S^„"= i/np" + aV, (1.10a) 
S^^" = -i/np" + alp". (1.10b) 

Pour conclure, la prescription envisagee reproduit parfaitement la decom- 
position de Haldane [Hal83a, Lev97]. 



2 Le « sur-reseau » antiferromagnetique 
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2.3 Cas d'un systeme frustre : le reseau triangulaire 

Le reseau triangulaire est pave avec le triangle equilateral represente 
dans Yalgebre reelle a divisions C par le triplet 

Ti = {l,w,w2}, (1.11) 
avec w la racine troisieme de I'unite : 

\22L 

w = e 3 . 



w 








^\J+1 


\ 






w — 





Figure 1.7: Representation complexe du triangle equilateral T . 



Le triangle equilateral est preserve par toutes combinaisons de I'operation 
de conjugaison et de la rotation w. Un « sur-reseau » etant choisi, les trois 
noeuds d'un «sur-noeud» quelconque (p) representes respectivement par les 
trois complexes 1, w et seront etiquetes respectivement par (po), (pi) et 
(^2)- L'antioperation antiferromagnetique n'etant rien d'autre que la rota- 
tion w, nous construisons le parametre d'ordre antiferromagnetique local rip" 
tel que 



V2 



(1.12a) 



et la representation triviale Ip" telle que 



(1.12b) 
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Soulignons que la representation complexe (1.12a) est definie a un champs de 
jauge locale 7p pres. A la limite classique les degres de liberie antiferromag- 
netiques satisfont aux contraintes suivantes : 



4 i2 



(1.13a) 
(1.13b) 



Le nombre total de degres de liberte (egal a six) est done preserve lors du 
changement des degres de libertes locaux et le parametre d'ordre antiferro- 
magnetique local rip"^ est effectivement norme a la limite continue. 

II s'avere alors judicieux d'introduire les combinaisons lineaires locales 



^Yp = ^ (Spi — ) ; 
a la limite classique nous avons 



(1.14a) 
(1.14b) 



(1.15a) 
(1.15b) 
(1.15c) 



Comme a la limite continue le couple l^rLxp",riYp"j n'est rien d'autre qu'un 
diedre orthonorme, introduisons un nouvel operateur n^p" tel que le triplet 



7^ ^p ) 



UYp" = 1+2^ - ^ lp\ 



("TLxp", TLyp", Tizp") forme un triedre orthonorme a la limite continue 



en injectant dans (1.16) les expressions (1.14) nous obtenons 
La representation complexe (1.12a) s'ecrit immediatement 
la representation pseudo-scalaire (1.16) se reecrit alors 



n " — 1 p-i7p 



Pa 



I e, 



a/37 



(1.16) 



(1.17) 



(1.18) 



(1.19) 



Consequence immediate : pour la theorie des champs ainsi construite le 
parametre d'ordre local vit dans SO (3), le groupe des rotations dans I'es- 
pace Euclidien R^. 



2 Le « sur-reseau » antiferromagnetique 
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Les spins internes d'une « sur-particule » (p) d'un reseau triangulaire an- 
tiferromagnetique sont donnes par les relations suivantes : 



C (X 


1 

V2 


V 




+ e" 




^1 ~ 


1 

V2 


V 




+ we~ 




C Ct 
•^2 ~ 


1 
V2 


V 


we+'>n/ 


+ we~ 





+ 0? I/, 



c 



(1.20a) 
(1.20b) 
(1.20c) 



En outre a chaque « sur-particule » (p) nous associons la chiralite Xp definie 
par [Siv95, Kaw98] 



(1.21) 



Un calcul rapide montre que la chiralite Xp est independante du champs de 
jauge locale 7p ; nous avons 



Xp — 17 ^Zp 1-pa- 



Cette relation (1.22) nous invite a introduire la chiralite nue 

Xp = ^Zp 



(1.22) 



(1.23) 



qui decrirait I'etat antiferromagnetique de la « sur-particule » : une « sur- 
particule » idealement antiferromagnetique (c'est-a-dire «a la Neel») aurait 
une chiralite nue nulle independamment du pas a du reseau et du nom- 
bre quantique s. La definition d'une chiralite nue Xp pour une « sur-particule » 
(p) appartenant a un reseau antiferromagnetique assouvi est immediate. 
Selon la relation (1.8b), les « sur-particules » de la chatne antiferromagne- 
tique ont une chiralite nue nulle. 



2.4 Cas d'un systeme assouvi : le reseau carre 

Le reseau carre est pave avec le 2-cube represente dans V algehre reelle 
a divisions C par le quadruplet 

Ci = {+l,+i,-l,-i}. (1.24) 



Le carre est preserve par toutes combinaisons de I'operation de conjugaison 
et de la rotation i. Un « sur-reseau » etant choisi, les quatre noeuds d'un «sur- 
noeud» quelconque (p) represente respectivement par les quatre complexes 
+1, -\-'\, —1 et — i seront etiquetes respectivement par (po), (pi), (P2) et (pa). 
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L'antioperation antiferromagnetique etant la rotation —1, nous construisons 
le parametre d'ordre antiferromagnetique local rip" tel que 



^ az/d/ = e" 



4 Z/TLp = [SpQ + ] ~ [%i + ] , 
la representation complexe dp" telle que 

[•Jpo ^ ' -^1 •->P2 ' •-'PS J 

et la representation triviale Ip" telle que 

4 ci Ip Spg -|- Sp^ -\- ~l~ . 
A la limite classique les contraintes sont : 



l-a^dp"tdp^ 



Tip^dp^ 



21 d "I 



TL "I 
' ^p '-p^ 



Vdp« + "^dp°"^dp^ 



(1.25a) 



(1.25b) 



(1.25c) 



(1.26a) 
(1.26b) 
(1.26c) 



Le nombre total de degres de liberte (egal a huit) est done preserve et le 
parametre d'ordre antiferromagnetique local Up" est effectivement norme a 
la limite continue. 

Selon la relation (1.26c) la chiralite nue Xp d'une « sur-particule » (p) du 
reseau carre est donnee par 



Xp 



dp"dp^ + Mp^Mp^ 



;i.27) 
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Aussi si la « sur-particule » (p) est idealement antiferromagnetique {chiral- 
ite nue nulle), alors dp" est un diedre orthogonal et le couple (dp", tip") 
un triedre orthogonal a la limite continue. Par consequent pour la theorie 
des champs ainsi construite le parametre d'ordre antiferromagnetique local 
pertinent vit dans S0(3), comme pour le reseau triangulaire. 

Les spins internes d'une « sur-particule » (p) d'un reseau carre antiferro- 
magnetique sont donnes par les relations suivantes : 



c ct 

^0 ~ 


tl . " 


+ 


75 «^ 




+ 


e~ 


>"fd 

U.p 


+ a' Ip", 


(1.28a) 


c a 

^1 ~ 


— i/TLp" 




75 «^ 






i e~ 


U-p 


+ anp", 


(1.28b) 


C CK 


I^TLp" 




75 «^ 


" e+'>d/ 


+ 


e~ 


>^dp"^ 




(1.28c) 


c en 


— i/Tlp° 


+ 


75 «^ 


ie+'>dp" 




i e~ 






(1.28d) 
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1 Le modele a non-lineaire 

Le modele a non-lineaire, ou la forme harmonique pour les mathematiciens, 
est un outil incontournable en physique theorique et un domaine de recherche 
en mathematique [Ura93, TN89, FraQl]. En physique de la matiere conden- 
see, le modele a non-lineaire decrit des materiaux antiferromagnetiques isotropes 
[FraQl, BP75, Tri79, CHN88, HalBSb]. 

1.1 Preliminaires 

Dans le cas present, le modele a non-lineaire correspondra a la limite 
continue du modele de Heisenberg decrivant un systeme de spins classiques 
distribues sur un reseau regulier bidimensionnel et interagissant entre plus 
proche voisin a priori. Le parametre d'ordre local [TK76] d'un tel systeme 
peut-etre represente par un vecteur unitaire fi generalement interprets dans 
la litterature comme une magnetisation locale [Nak90] : dans le Chapitre 
precedent, je montre que la nature du parametre d'ordre n est plus complexe. 




Figure 2.1: Le parametre d'ordre fi : coordonnees spheriques. 
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Figure 2.2: Le parametre d'ordre n : coordonnees stereographiques. 

Le Hamiltonien magnetique J^mag d'une membrane magnetique Ai sera 
done decrite par la fonctionnelle [DVGSB95, VGDSB95, BD98] 

^mag = J j j ^ g'^ K^diu"^ d.n'' , (2.1) 

M 

ou I'integrale est prise sur toute la membrane J^. La constante J designe le 
couplage spin-spin isotrope entre plus proches voisins, c'est-a-dire le terme 
d'echange electronique [Lev97] entre deux sites voisins du reseau initial. Les 
tenseurs metriques gij et hap decrivent respectivement la membrane magne- 
tique Ai et la variete du parametre d'ordre n, autrement dit la sphere unite S'^. 

Le parametre d'ordre n, qui doit satisfaire la contrainte = 1, admet 
deux descriptions naturelles : une description geometrique et une description 
algebrique. La description algebrique est implicitement utilisee dans le Hamil- 
tonien (2.1), toutefois un systeme de coordonnees adequate reste a preciser. 
Le parametre d'ordre n est communement represente par ses coordonnees 
spheriques (6,$) (Figure 2.1) : 

n=[sinO cos$,sinO sin$,cos0]. (2-2) 

Le parametre d'ordre fi est ainsi represente par un point sur la sphere unite S'^ ; 
par consequent la metrique h de la variete cible s'ecrit 

/i = de^de + sin^e d$®d$. (2.3) 
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Le Hamiltonien magnetique (2.1) devient alors 




(2.4) 



M 



Sous cette forme, qui incorpore la contrainte du parametre d'ordre, le carac- 
tere non-lineaire de la fonctionnelle (2.1) apparait ostensiblement. 

La description algebrique consiste a projeter le parametre d'ordre n sur 
la sphere complexe C U {oo} via une projection stereographique [Coh67, 
Nak90]. Rappelons au passage que la sphere complexe (ou de Riemann) 
s'obtient par compactification du plan complexe C : I'ensemble des points 
situes a I'infini est identifie a un point unique, I'infini. Compte tenu des 
conventions precedemment employees, il est commode de projeter stereo- 
graphiquement la sphere unite ^S"^ par rapport au pole sud. Les coordon- 
nees stereographiques (U, V) associees au parametre d'ordre n correspondent 
alors aux coordonnees cartesiennes du point intersection du plan equatorial 
et de la demi-droite issue du pole sud passant par le point de coordonnees 
spheriques (0, $), representant du parametre d'ordre n sur la sphere unite S"^ 
(Figure 2.2). Notons que le pole nord (© = 0) se projette sur I'origine, les 
points de I'equateur (© = 7r/2) sur eux-meme, et le pole sud (© = tt) sur 
le point infini. Des notions elementaires de geometric analytique conduisent 
aux relations de passages suivantes : 



U = tan i© cos 
V = tani© sin$. 



(2.5a) 
(2.5b) 



Au parametre d'ordre n est alors associe la variable complexe 



W = U + \V 



(2.6a) 



qui vit effectivement sur la sphere complexe C U {oo} ; nous lisons 



=tani©e'*. 



(2.6b) 



Ici, le Hamiltonien magnetique (2.1) s'ecrit immediatement 




(2.7) 



en incorporant dans (2.4) I'egalite utile 



sine diW^W [diO + isin© 9^$] . 



(2.8) 
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1.2 Double obstruction topologique : 
classification de Toulouse-Kleman 

Considerons maintenant la membrane M. comme rien d'autre qu'une var- 
iete orientable de dimension 2 fermee ou ouverte qui, si necessaire, peut- 
etre compactifiee en identifiant chacun de ses bords, si cela a un sens, a 
un point unique ; nous noterons M la membrane Ai ainsi compactifiee. En 
d'autres termes, la membrane est momentanement envisagee comme un objet 
topologique : seules nous interessent ici les proprietes qui ne changent pas lors 
de deformations continues. 

Chaque membrane compactifiee M est ainsi classee suivant son genre 
topologique g G N : deux membranes compactifiees M de meme genre topologique 
g sont deformables continument de I'une a I'autre sans introduire aucune 
dechirure ni trou [Nak90]. Pour illustration : une sphere est de genre topologique 
0, un tore de genre topologique 1, deux tores sondes de genre topologique 2, 
etc. Autrement dit, I'invariant topologique g verrouille le nombre de trous de 
la membrane ; nous noterons Tg la classe d'equivalence de genre topologique 
g. La topologie combinatoire fournit une relation simple entre la caracteris- 
tique d'Euler et le genre topologique [Nak90] : 



tandis que le theoreme de Gauss-Bonnet affirme [Str61, DFN92, Nak90] : 



avec K la courhure locale de Gauss (2.28). 

Une approche similaire existe pour les configurations de spins. En efi^et, 
les configurations continues du parametre d'ordre n ne sont autres que les 
fonctions continues qui envoient la membrane compactifiee M (la variete sup- 
port) sur la sphere unite (la variete cible). Or I'ensemble de telles fonc- 
tions forme le groupe de cohomotopie tt^ (Tg) , avec g le genre topologique 
de M : les fonctions continues sont classees ici selon leur classe d'homotopie 
[Spa49, HuoQ, Xak90]. Consequence immediate : si le groupe de cohomotopie 
TT^ (Tg) ne se reduit pas au groupe trivial, nous nous attendons a rencon- 
trer une physique exotique. En elfet, certains concepts generalement tres 
feconds en physique sont ici infructueux : les perturbations « continues » ou 
« plastiques » ne permettent pas de sauter d'une classe de configuration a une 
autre. Aussi, la classifi,cation des configurations de spins suivant leur classe 
d'homotopie s'avere fort legitime, voire essentielle [BP75, Bog76, TK76]. Les 
configurations de spins homotopes a la configuration de spins dont la fonc- 
tion continue associee n envoie la membrane compactifiee M en un point 



x(Tg) = 2(l-g); 



(2.9) 




(2.10) 



M 
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unique de la sphere unite ^S"^ {e.g. au pole nord) correspondent aux configu- 
rations triviales. En fait, I'homologie et la cohomologie nous apprennent que 
les groupes de cohomotopie tt^ (Tg) sont isomorphes au groupe des entiers 
relatifs Z. Nous avons effectivement 

(^2^ = 772 (5') ^ Z, (2.11) 
et, en invoquant la formule de Runneth, 

TT^ (r^) = TT^ X S') = tt' (s') ® 7i' (S') ^ Z; (2.12) 
done, par sommation connexe, 



TT^ (Tg) ^ Z Vg e N. (2.13) 

U invariant topologique qui verrouille les configurations de spins est Yindice 
de Pontrjagin Q qui s'ecrit 

Q = ^JJ sine ded^ (2.14) 

M 

dans le contexte present [Fel87, Nak90]. Contrairement aux topologues, les 
physiciens parle de charge topologique ; nous noterons tt^q la classe d' equiv- 
alence de charge topologique Q qui envoie la classe d'equivalence Tg sur la 
sphere unite 5*^. Get invariant topologique correspond an nombre de fois que 
la sphere unite S"^ (la variete cible) est recouverte entierement quand la mem- 
brane compactifiee M (la variete support) est parcourue [Fel87]. En adoptant 
le point de vue originelle, Yindice de Pontrjagin decompte les retournements 
du parametre d'ordre ii. Par consequent, les configurations triviales de spins 
seront caracterisees par une charge topologique nuUe ; par contraste, les con- 
figurations topologiques auront une charge topologique non-nuUe. 

Pour resumer, la membrane magnetique M. est caracterisee par deux in- 
variants topologiques : 

- son genre topologique g qui verrouille sa configuration spatiale ; 

- sa charge topologique Q qui verrouille sa configuration magnetique. 

II est evident qu'une modelisation pertinente de la membrane magnetique 
doit a la fois satisfaire et clarifier cette double obstruction topologique dont 
le mecanisme peut-etre schematise comme suit : 

, , compactification , , ~ 7r^(rg)=Z ^ 9 ^ _\ 

M > M e Tg — > n e TT^Q. (2.15) 

genre topologique g charge topologique Q 
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1.3 Energie magnetique : 

decomposition de Bogomol'nyi 

Comme des considerations de symetrie permettent de modeliser de nom- 
breux systemes physiques a ordre symetrique et d'elaborer des theories physiques 
[Siv95], des considerations de topologie doivent permettre de decrire des sys- 
temes physiques a ordre topologique [TK76, Mer79, Siv95]. Notons en re- 
vanche que, si I'origine des invariants de symetrie est acquise, I'emergence 
invariants topologiques en physique demeure une enigme [Ryd85]. II est 
done remarquable que nous puissions retrouver le verrouillage topologique in- 
duit par la distribution de spins a partir du Hamiltonien magnetique (2.1) en 
appliquant une methode suggeree independamment par Belavin et Polyakov 
[BP75] d'une part, et par Bogomornyi [Bog76] d'autre part. 

Pour commencer, introduisons le tenseur T^" tel que 

T,- = l=[d,n--sE;F'^,drn^] £ = ±1. (2.16) 

Les tenseurs Eij et Fa/s sont les tenseurs unitaires parfaitement antisymetriques 
d'ordre deux associes respectivement a la membrane magnetique et a la 
variete du parametre d'ordre n ; nous lisons 

Eij = y/g Sij, Vheaf), (2.17) 

avec Cab le pseudotenseur unitaire parfaitement antisymetrique. En choisis- 
sant pour la variete cible la metrique spherique (2.3) et en reduisant locale- 
ment la metrique g a une metrique orthogonale, il saute aux yeux que le 
scalaire g'^^hapT^'^Tj^ est positif. Or, puisque Tj" est un tenseur, cette pro- 
priety reste vraie pour tout autre choix de systeme de coordonnees. Nous en 
deduisons que 

g'^h^fsT.'^T/^O. (2.18) 



En outre, le calcul explicite de cette trace donne 



g'^KpT.^T/ = g'^Kpdiu'^djn^ - 2e^J, (2.19a) 

V 9 

ou J est le Jacobien de la transformation locale {x^) — > {rf) : 

J ^\e,xe''drn'^dsn\ (2.19b) 
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Des relations (2.18) et (2.19a) nous tirons la precieuse inegalite 



\/9 



(2.20) 



L'integration membre a membre de I'inegalite (2.20) sur toute la membrane 
compactifiee M en adoptant la metrique spherique (2.3) pour la variete cible 
nous conduit aisement au point remarquable de la methode : 



JJ^dQ g'^hapd,nf'djn^ ^2e JJVhdQ J, 

M M 



>2 



jj sin e de d$ 

M 

^SttIQI, 



(2.21a) 

(2.21b) 
(2.21c) 



avec Q Yindice de Pontrjagin (2.14). Ayant encore en memoire I'expression 
du Hamiltonien magnetique (2.1), nous ecrivons sans attendre 



'^maq ^ SttJ IQI 



(2.22) 



Plus prosaiquement, I'energie magnetique des configurations de spins ap- 
partenant a une meme classe d' equivalence TTgQ est minoree par Venergie 
topologique 



Eq = 87rJ|Q| . 



(2.23) 



Qui plus est, par construction, cette energie topologique Eq n'est atteinte que 
par les configurations de spins pour lesquelles le tenseur Tj" se confond avec 
le tenseur nul : c'est-a-dire, compte tenue de la definition (2.16) du tenseur 
pour lesquelles le parametre d'ordre fi verifie les equations 



(2.24) 



Si la metrique de la variete cible est la metrique spherique (2.3), les equations 
auto-duales (2.24) se lisent 



diO = e ^/g sin 6 eird^^ ; 



(2.25) 



tandis que leur versions algebriques s'ecrivent 



(2.26) 
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Les configurations topologiques auto-duales seront done les configurations topologiques 
qui minimisent V energie magnetique (2.1). Enfin, Yenergie topologique Eq 
s'interprete immediatement comme I'energie minimale necessaire pour re- 
tourner |Q| fois le parametre d'ordre n, et le facteur SttJ comme le quantum 
d'energie minimale de retournement. Avant de conclure, notons que I'injec- 
tion de I'egalite (2.19) dans la fonctionnelle originelle (2.1) donne une ecriture 
tres pertinente du Hamiltonien magnetique [Fel87] : 

J^mag = J J J ^ dfl g'^ h^^p T,"T/ + 87rJ|Q| . (2.27) 

M 

Pour resumer, la decomposition de Bogomol'nyi montre que le Hamiltonien 
magnetique (2.1) live la classification topologique des configurations magne- 
tiques : 

- I'energie magnetique des configurations de spins appartenant a la meme 
classe d'homotopie est minoree par I'energie minimale necessaire pour 
retourner le parametre d'ordre n ; 

- les configurations topologiques qui minimisent I'energie magnetique 
verifient les equations auto-duales (equations differentielles du premier 
ordre) . 

Un mecanisme comparable existe pour la classification topologique des con- 
figurations spatiales. 



1.4 Energie de courbure : 

Hamiltonien de Helfrich-Willmore 

Pour commencer, rappelons qu'une surface bidimensionnelle est parfaite- 
ment decrite localement par ses deux courbures principales Ai et A2. Leur 
produit 

K = A1A2 (2.28) 
est appele la courbure locale de Gauss de la surface, leur moyenne arithmetique 

H = ^[X, + A2] (2.29) 
la courbure moyenne locale de la surface. 
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Figure 2.3: Le tore de Clifford est le tore a symetrie axiale dont le rapport 
du rayon axiale r sur le rayon de revolution R est egal a l/V^. Selon une 
conjecture due a Willmore [Wil82], \a fonctionnelle de Willmore est minimisee 
soit par le tore de Clifford soit par un de ses transformes conformes. 



Le Hamiltonien gouvernant la dispersion des courhures principales au sens 
de la theorie des probabilites [Pin86, Wil93] 



M 



Ai — A2 



n 2 



(2.30) 



apparait comme un candidat legitime pour decrire les deformations continues 
et localisees subies par une membrane souple compactifiee M spontanement 
a courbure constante. La constante ko sera interpretee comme la rigidite de 
courbure nue. La relation immediate 



Ai — A2 


2 


"A1 + A2' 


2 




2 



— A1A2 



(2.31a) 
(2.31b) 



permet d'ecrire le Hamiltonien elastique naff (2.30) sous la forme plus perti- 
nente : 



(2.32) 



M 



M 
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D'apres le theoreme de Gauss-Bonnet (2.10), I'integrale de droite ne depend 
que du genre topologique g de la membrane souple compactifiee M. L'integrale 
de gauche n'est autre que la fonctionnelle de Willmore [Wil82, Wil93] : pour 
des surfaces a topologie spherique (g = 0), elle est minimale (et nulle) pour la 
membrane plane; pour des surfaces a topologie torique (g = 1), elle est mi- 
noree par An et serait minimale {conjecture due a Willmore [Wil82]) soit pour 
le tore de Clifford (FIGURE 2.3) soit pour un de ses transformes conformes. 
Ainsi le Hamiltonien elastique nai'f (2.30) ne convient-il plus pour decrire 
des deformations continues et localisees endurees par une membrane sou- 
ple compactifiee M de topologie non spherique (g 7^ 0). En introduisant une 
courbure moyenne spontanee non nulle Hq, Helfrich [Hel73, Fou92] obtient 
un Hamiltonien elastique plus adequate : 

je,i = \Kjj^dn [H-Hof-lk,JJ^dQK, (2.33) 

M M 

ou la rigidite de courbure nue ko a ete renormalisee. La constante kc designe 
la rigidite de courbure, kg la rigidite topologique. Dorenavant, par Hamiltonien 
elastique J^^i nous entendrons le Hamiltonien heuristique (2.33). En vertu du 
theoreme de Gauss-Bonnet (2.10) et de la relation (2.9), le Hamiltonien de 
Helfrich-Willmore (2.33) se lit 

^ei = Ik^JJ^dn [H - Hof + 27rA;g (g - 1) . (2.34) 

M 

Cette ecriture du Hamiltonien elastique M'^ doit etre comparee avec I'ecrit- 
ure (2.27) du Hamiltonien magnetique Jifmag, et reciproquement. En regard 
de quoi, la courbure spontanee Hq qui doit decrire les proprietes physiques de 
la membrane souple A4 correspondra dans la suite du Chapitre a la courbure 
moyenne d'une surface choisie a priori arbitrairement comme etant effec- 
tivement la forme spontanee de la membrane compactifiee M. Dans ce cas 
restrictif, le Hamiltonien magnetique (2.34) se contente de decrire des defor- 
mations de faible amplitude, et les courbures renormalisees kc et kg dependent 
clairement de ce choix arbitraire. 

Pour resumer, nous pouvons construire un Hamiltonien elastique qui live 
la classification topologique des configurations spatiales : 

- V energie elastique est minoree par I'energie minimale requise pour creer 
un trou topologique ; 

- les formes topologiques qui minimisent I'energie elastique correspondent 
aux formes spontanees de la membrane compactifiee M. 
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1.5 Frustration geometrique 

Dans les parties anterieures nous montrons qu'une membrane magnetique 
M. se caracterise par une double classification topologique, chacune etant 
levee par un Hamiltonien : 

- le Hamiltonien elastique Jifei engendre des formes homeomorphes a des 
formes spontanees pour lesquelles V energie de courbure E^i est saturee ; 

- le Hamiltonien magnetique M'mag genere des confi,gurations de spins 
homotopes a des configurations topologiques auto-duales pour lesquelles 
V energie magnetique E^q^ est saturee. 

Ainsi a la double obstruction topologique (2.15) correspond-t-il la double 
levee topologique suivante : 

^ ^ ^ ^ ^ ^ g ^ o _ ^2.35) 

forme spontanee configuration auto-duale 

Generalement, le Hamiltonien total 3^ d'une membrane magnetique M. 
est defini comme la somme de trois Hamiltoniens : 

= t^^'mag ~l~ '^^el ~l~ '^^m—eh (2.36) 

ou J^mag, J^ei Gt .^m-ei repiesentent, respectivement, le Hamiltonien mag- 
netique, elastique et magneto-elastique. Or, pour un systeme quasi-unidimen- 
sionnel, I'addition d'un terme de couplage spin-elasticite revient a renormaliser 
la constante J de couplage spin-spin [CT80]. Nous definirons done le Hamilto- 
nien total de la membrane magnetique compactifiee M comme la somme 
du Hamiltonien magnetique (2.1) et du Hamiltonien elastique (2.33) ; des 
expressions (2.27) et (2.34) nous tirons la formule 

^ = kcffVgdn In g'^KpT^'^T/ + \[H - Ho[ 



M 



+ 87rJ|Q| +27rA;g(g- 1) , (2.37) 



ou nous avons introduit la constante relative de couplage n = J /kc- L'ecriture 
(2.37) du Hamiltonien total place pertinemment sur un meme pied les 
interactions elastique et magnetique : la metrique gij et le parametre d'ordre 
n interagissent. Par configuration geometrique nous designerons un couple 
(5fjj,n) qui minimise le Hamiltonien total J^, c'est-a-dire un couple {gij^rv) 
solution des equations d'Euler-Lagrange deduites du Hamiltonien Ainsi la 
double levee topologique (2.35) se combine-t-elle en une levee double : 

M > {gij,n) e Tg X Tr^Q. (2.38) 

configuration gfeomfetrique 
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A priori les configurations geometriques (gij, n) ne saturent pas le Hamil- 
tonien total M' : les configurations geometriques dont I'energie totale E est 
superieure a Venergie minimum minimorum topologique 

E = 87rJ|Q| +27rA;g(g- 1) (2.39) 

sont dites frustrees. Cette frustration geometrique apparait clairement comme 
la competition entre les deux ordres topologiques, et I'exces d'energie qui la 
caracterise comme une energie interaction topologique. 

Pour resumer, le Hamiltonien total Jif combine la double classification 
topologique de la membrane magnetique souple compactifiee M : 

- I'energie des configurations geometriques appartenant a la meme classe 
d'homotopie est minoree par Venergie minimum, minimorum topologique 
E necessaire pour retourner le parametre d'ordre n et pour former les 
trous topologiques ; 

- une frustration geometrique emerge a priori de la competition entre 
les deux ordres topologiques, les configurations geometriques (qui veri- 
fient les equations d'Euler-Lagrange deduites de J^) ne minimisant pas 
necessairement a priori Venergie magnetique et/ou Venergie elastique. 

II est evident que des illustrations simples doivent permettre de clarifier le 
mecanisme de cette frustration geometrique. 
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2 Membranes de topologie spherique 

Dans cette partie, nous envisageons des membranes magnetiques souples 
compactifiees de genre topologique spherique et de forme spontanee simple : des 
brisures de symetrie successives revelent une physique de plus en plus exotique. 

Je met en evidence le mecanisme de la frustration geometrique en gelant une 
configuration magnetique de la membrane souple rigidifiee. 



2.1 La sphere : une forte degenerescence 



Pour une membrane souple rigidifise {kc = oo), le Hamiltonien total 
se reduit au Hamiltonien magnetique J^mag '■ les configurations geometriques 
se confondent alors avec les configurations topologiques auto-duales. 

Une representation naturelle de la sphere rigide, en coordonnees spheriques 
(p, e,(p), est 



(2.40) 



avec r le rayon de la sphere rigide. La metrique de la membrane magnetique 
s'ecrit alors 



g = r 



de^de + 8111^ e df^od^? 



ainsi g^'^ — g'^^ = et nous avons 

/^V^ = sin^, 



sin^ 



(2.41) 



(2.42) 



Aussi I'expression (2.4) du Hamiltonien magnetique s'ecrit 

+ 



-^mag = J j dO j dip 

— TT 



. _ ^2 sin^ 6 ^2 
sm 9Ql + — — $^ 



sin^ 



+sin^ sin^e $2 



sin 9 



(2.43) 



ou les indices representent les derives partielles. Notons que la densite Hamil- 
tonienne magnetique de la sphere rigide magnetique est independante du 
rayon r de la sphere. De la fonctionnelle (2.43) nous tirons aisement les equa- 
tions auto-duales de la sphere rigide magnetique : 



sin^ Qe — £ sin© 

O,^ = —e sm.9 sin© ^q. 



(2.44a) 
(2.44b) 
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Leur versions algebriques prennent la forme 



(2.45a) 
(2.45b) 



Les nouvelles coordonnees (m, v) correspondent aux coordonnees stereographiques 
par rapport au pole sud de la sphere rigide projetee sur la sphere unite ; 
nous lisons 



u — tan ^9 cos ip, 
V = tan sin (p. 



(2.46a) 
(2.46b) 



Reconnaissant le critere de Cauchy-Riemann , nous affirmons que W est une 
fonction meromorphe [WWo2] de la variable complexe w = u + e\v, le signe e 
correspondant au signe de la charge topologique Q. Par consequent, il existe 
une infinite de configuration topologique auto-duale sur la sphere magnetique 
rigide ; la configuration topologique auto-duale W {w) — w est communement 
appelee la « configuration herisson » [Fra91, Tho98]. 

Par la suite, la recherche de configurations symetriques permet de main- 
tenir les calculs abordables. Dans le cas present, il s'agit de choisir arbitraire- 
ment un des axes de rotation de la sphere spontanee puis de rechercher les 
configurations symetriques suivant cet axe. 

Compte tenu des conventions deja utilisees et en introduisant la nouvelle 
coordonnee ( = 7r/2 — 6 afin de souligner la symetrie par rapport au plan 
equatorial, les configurations symetriques sont ici telles que 



$^ = 



avec comme conditions aux poles 

©(C=-f) = et e(C-+f)-0[7r] 
Le Hamiltonien magnetique (2.43) s'ecrit alors 



_i -TT 



cos C ©J + 



9 sm 6 ^9 



cosC 



(2.47) 
(2.48) 

(2.49) 



Les equations d'Euler-Lagrange des configurations symetriques qui extremisent 
le Hamiltonien magnetique {SM^mag — 0) sur la sphere rigide magnetique se 
deduisent de (2.49) sans diflficulte : 



^ (pep 

cos C d^[cos C ©c] — ® ® ^¥ 



(2.50a) 
(2.50b) 
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L'integration de I'equation (2.50a) est immediate : 

= ciif q<^ e (2.5ia) 

par consequent (2.50b) devient 

cos C df [cos C 6^] = sin 6 cos 6. (2.51b) 

Cette nouvelle equation (2.51b) s'integre sans effort en multipliant ses deux 
membres par ; nous lisons 

cos^ C ^ q J sin^ © + ^ me[0,+oo[. (2.o2a) 

Or les conditions aux poles (2.48) imposent 

m = 0, (2.52b) 

aussi nous aurions du lire 

cosC ©c — ^ q^ ^i^®- (2.52c) 

Consequence immediate : les configurations auto-duales recouvrent toutes 
les configurations magnetiques symetriques puisque (2.51a) et (2.52c) veri- 
fient sans equivoque I'equation auto-duale (2.44a). Des calculs analytiques 
fastidieux montrent que Vecriture algebrique des configurations magnetiques 
auto-duales est de la forme 

W{w) = w'^^. (2.53) 

Pour £q^ = 1, nous retrouvons la « configuration herisson» 

e = f - C = ^. (2.54) 

En outre, des calculs ennuyeux de geometric differentielle montrent que les 
deformations symetriques de la sphere souple sont d'ordre trois au moins : 
les configurations magnetiques symetriques ne deformerons done pas une 
sphere souple suffisamment rigide (k ~ 1). Autrement dit, les configurations 
geometriques symetriques de la sphere souple magnetique ne sont pas frus- 
trees . 
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2.2 Le cylindre infini : des solitons assouvis 



Considerons maintenant une membrane magnetique souple compactifiee 
M toujours de genre topologique spherique mais dont la forme spontanee 
n'admet qu'un seul axe de revolution : une section de cylindre droit. 

En coordonnees cylindriques {p, (p, z) , une telle membrane rigidifiee est 
parfaitement decrite par 



p = r, 
-Az ^ +Az, 



(2. .55a) 
(2.55b) 



ou les constantes reelles positives r et Az sont respectivement le rayon du 
cylindre et la demi-longueur de la section. La metrique de la membrane mag- 
netique est alors de la forme 



g — r^d(/7(8)d(/7 + dz<^dz; 
aussi g'^^ — g^'^ = et nous avons 

D'emblee nous etudions les configurations symetriques telles que 

0; 

dans ces conditions le Hamiltonien magnetique (2.4) s'ecrit 



+A2 +7r 



-Az -TT 



^2 sin^ ^2 



(2.56) 



(2.57) 



(2.58) 



(2.59) 



Une expression sans dimension de (2.59) s'obtient en effectuant un change- 
ment d'echelle suivant I'axe des z : 



^ =7 



AC -TT 



0^ + sin^ $2 



avec 



C = z/r et = Az/r. 



(2.60) 



(2.61) 



Enfin, la membrane magnetique est compactifiee en imposant une densite 
Hamiltonienne nulle aux bords : 



0C (C = ±AC) = 0, 
0(C = -AC)=O et 0(C = +AC) = O[7r]. 



(2.62a) 
(2.62b) 
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Les equations auto-duales symetriques d'une section du cylindre magnetique 
rigide se deduisent facilement de la fonctionnelle (2.60) : 

= £ sine (2.63) 

En outre les relations (2.58) imposent le dedoublement 

= q^, (2.64a) 
= £ sine, (2.64b) 

avec q<^ G Z. Enfin la solution de I'equation (2.64b) est le soliton 

e (C) = 2 arctan [exp {eq^Q] . (2.65a) 

Par consequent les conditions aux bords (2.62) ne sont satisfaites que si 

AC = +00, (2.65b) 

les solitons (2.65a) recouvrant alors jq^,] fois la sphere unite S'^. En d'autres 
termes, seul le cylindre de longueur infinie admet des configurations magne- 
tiques auto-duales symetriques. Avant d'approfondir le cas des cylindres de 
longueur finie, achevons celui du cylindre infini. 

Puisque, compte tenu du changement d'echelle (2.61) et de I'expression 
des solitons magnetiques (2.65), le rayon r et la coordonnee z apparais- 
sent respectivement comme le parametre de forme pertinent et la coordonnee 
curviligne pertinente, le rayon r de la membrane souple rigidifiee doit etre 
libere suivant I'axe des z : 

r(z)=ro [1 + A{z)], (2.66) 

ou To represente le rayon spontane et la fonction A la deformation continue 
locale. L'orthogonalite de la metrique rigidifiee est conservee tandis que les 
formules (2.57) deviennent 



1 + r^^l r 
f^^r9 = ' / . 9''^^9 = I • (2-67) 



Considerant des deformations continues de faible amplitude, nous develop- 
pons la densite du Hamiltonien elastique (2.34) jusqu'a I'ordre deux en A, 
K et A^z : 

J^i^lTTkc JdCA'', (2.68) 

-AC 

avec des notations naturelles et en omettant le terme topologique. 
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Si une des configurations magnetiques de la membrane souple M rigid- 
ifiee est gelee, le Hamiltonien magnetique J^mag de la membrane souple M 
devient alors une fonctionnelle de la seule deformation A. Dans ce contexte 
les deformations A qui minimisent le Hamiltonien total ^ permettent de se 
faire une idee des deformations effectivement induites par la distribution de 
spins et subies par la membrane souple compactifiee M. 

Pour la metrique souple (2.67), un developpement limite de la densite du 
Hamiltonien magnetique (2.4) jusqu'a I'ordre deux en A et A^ donne : 



+AC +7r 

-AC -TT 



+ sin^ e $J 



(-A + 1A9 



- sin^ e $2 



+ A^ sin^ e 



, (2.69) 



pour une configuration magnetique symetrique gelee quelconque. Des lors les 
deformations A induites par les solitons magnetiques (2.65) geles et subies 
par le cylindre souple magnetique infini doivent extremiser le Hamiltonien 
total 



jr = 27rJqjy'dC 



2 + A^ 



sin^ e + 



1 



16«;q2 



A^ 



(2.70) 



somme du Hamiltonien elastique (2.68) et du Hamiltonien magnetique (2.69) 
pour ces configurations magnetiques. Or il apparait clairement que seule la 
deformation triviale A = minimise cette fonctionnelle (2.70) : les solitons 
magnetiques (2.65) ne deforment pas le cylindre souple magnetique infini. 
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2.3 Sections de cylindre : deformation solitonique 

Dans la section precedente, nous montrons qu'il n'existe pas de configu- 
ration magnetique symetrique auto-duale sur les sections finies de cylindre 
magnetique rigide : ce n'est pas pour autant qu'il n'existe pas de configuration 
topologique symetrique. 

De telles configurations doivent satisfaire les conditions aux bords (2.62) et 
minimiser le Hamiltonien magnetique (2.60), c'est-a-dire verifier les equations 
d 'Euler-Lagrange 

%^ = 0, (2.71a) 

e^^ = sin e COS e (2.71b) 

L'equation (2.71a) s'integrant sans efi'ort, nous ecrivons 

= q^p, (2.72a) 
e^c = qj, sine cose, (2.72b) 

avec e Z le nombre de recouvrement de la sphere unite S"^ autour de 
I'axe de revolution. Notons an passage que ce systeme dedouble (2.72) re- 
produit le dedoublement auto-dual (2.64). De surcroit l'equation (2.72b), au 
changement d'echelle g = pres, n'est rien d'autre que Vequation simple 
de sinus-Gordon (sg) ^ ; la solution de l'equation (2.72b) est done le soliton 
periodique 

e (C) = ea (q^C I ^) - f 5paire,qc, (2.73a) 

avec q^ e Z le nombre de recouvrement de la sphere unite le long de I'axe 
de revolution. Les conditions aux bords (2.62b) ne sont satisfaites que si le 
parametre m verifie 



m = ^ ^Cj , (2.73b) 

les solitons periodiques (2.73a) ayant alors la charge topologique 

Q = q^qc- (2.73c) 

Observons pour finir que les resultats du cylindre rigide infini se retrouvent 
aisement en faisant tendre vers I'infini, tandis que le rayon r et la coor- 
donnee z gardent leur pertinence. 

^Le lecteur est invite a consulter 1' Annexe A pour se familiariser (rapidement) avec 
Vequation de sinus- Gordon et les notations utilisees. 

Neanmoins, pour le lecteur averti : a{- \ fh) correspond a la solution croissante de Vequation 
simple de sinus-Gordon (sg) localisee a I'origine et caracterisee par le parametre rh, la 
quasi quart-periode de sinus-Gordon Ka verifiant Ka {fh) = K(l + m) avec K I'integrale 
elliptique de premiere espece. 
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Figure 2.4: Deformation symetrique d'une sec- 
tion de cylindre magnetique souple en presence 
d'un 27r-soliton magnetique (2.73) : I'unite de 
longueur est le rayon r de la section rigidifiee 
(representee par la grille exterieure), la demi- 
longueur relative = Az/r est egale a 5, la 
constante relative de couplage k = J/kc a 1/2, 
enfin la deformation locale A est magnifiee par un 
facteur 10. 

Le soliton frustre se libere en retrecissant globale- 
ment la membrane magnetique ; les deux renfle- 
ments correspondent aux deux vagues du soliton : 
son energie magnetique Emag s'y concentrant in- 
dependamment de la geometric de la membrane, 
il refoule I'energie magnetique accumulee en ses 
branches. 
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Afin d'etudier la frustration geometrique des solitons magnetiques (2.73), 
introduisons S'mag le rapport de leur energie magnetique E^^g sur leur energie 
topologique Eq (2.23) : 



E, 

jp — 

^mag — j-, 
ill, 



mag 



(2.74) 



Le calcul explicite de S'mag donne 



SI 



mag 



E(l + m) + im K(l 



m 



E 



1 



1 + m 



^ 1 + m 



1 + m 



(2.75a) 
(2.75b) 



ou les fonctions K et E designent respectivement I'integrale elliptique de 
premiere et de seconde espece [WW52, Law89]. L'energie relative Smag des 
solitons magnetiques (2.73) ne depend ainsi que du parametre m. Qui plus est, 
^mag croit strictement de 1 vers oo lorsque m varie de vers oo. Le parametre 
m se revele done comme le parametre pertinent de la frustration geometrique : 
plus m est grand plus les solitons magnetiques (2.73) sont frustres, le cas 
m = correspondant a un systeme assouvi. Retrospectivement la relation 
(2.73b) permet de mieux saisir le mecanisme de frustration des sections finies 
de cylindre magnetique rigide : les solitons magnetiques (2.73) tendent a se 
retourner autour de I'axe de revolution (|q^| grand) et a s'etendre le long 
de I'axe de revolution (Iq,;^! petit et grand). Consequence immediate : les 
solitons magnetiques (2.73) sont frustres parce que confines {A( < oo). 

Maintenant les deformations A subies par les sections de cylindre souples 
magnetiques et induites par les solitons magnetiques frustres (2.73) geles 
doivent extremiser le Hamiltonien total 



+AC 



27rJqJ JdC 

-AC 



(2 + m) - m -A + ^A^ 



-{1 + m) 2 + A^ sn^ (q^C | 1 + m) + 1 + 



16fi;q2 



lA^ 



, (2.76) 



qui doit etre lu comme la version frustree du Hamiltonien total (2.70). Les 
equations d'Euler-Lagrange de cet Hamiltonien (2.76) prennent la forme 



Agg + (1 + m)A — mB sn^ {g \ m) A — ^Jm j, 



(2.77a) 



2 Membranes de topologie spherique 



41 



ou nous avons pose 



m 
A 

B 



1 + m, 
2 



mq 



1 



-1 



(2.77b) 
(2.77c) 

(2.77d) 

(2.77e) 



Ici nous avons egalement introduit la nouvelle variable g = q^(. L'equation 
difFerentielle lineaire inhomogene du second ordre (2.77) n'est rien d'autre 
y equation inhomogene de Lame (il) avec un terme inhomogene constant ^ ; une 
solution raisonnable de l'equation (2.77) est done la deformation solitonique 



A(C) = L(q^C|l + m;A,5,j) 



(2.78) 



La frustration geometrique apparait ostensiblement dans l'equation d'Euler- 
Lagrange (2.77) : A = n'est pas solution. Concentrons-nous a present sur 
le mecanisme de deformation. D'apres les relations (2.61), (2.73b) et (2.75) 
nous Savons que Venergie magnetique Emag des solitons (2.73) decroit en 
fonction du rayon r : les solitons tendent a retrecir radialement la section de 
cylindre souple magnetique. D'un autre cote le Hamiltonien elastique (2.68) 
essaie de maintenir sa forme spontanea. Par consequent la competition entre 
Venergie magnetique et Venergie elastique se traduit par une reduction du 
rayon r. Or Venergie du soliton est essentiellement localisee sur sa vague (0 
proche de |) selon (2.60) : la deformation est ainsi moins importante sur 
la vague du soliton puisque I'energie elastique est isotrope. Pour interpreter 
les consequences geometriques de cette liberation geometrique definissons la 
dilatation relative du cylindre Xc par 



To 



1 + A. 



(2.79) 



^Le lecteur est encourage a lire I'Annexe B qui aborde Vequation inhomogene de 
Lame (il) en se bornant au contexte rencontre. 

Toutefois, pour le lecteur averti : L (• | m;A,B,j) designera la fonction (bornee) asso- 
cicc a la solution minimalc dc l'equation aux differences lineaire du second ordre dcduite 
de Vequation inhomogene de Lame (il) suivant la methode communement utilisee pour 
construire les fonctions dites de Lame (solutions de Vequation homogene de Lame). 
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Figure 2.5: La dilatation relative Ac de sections de cylindre magnetiques 
souples en presence d'un vr-soliton magnetique (2.73) : la dilatation rela- 
tive Ac definie par (2.79) est representee en fonction de (/A( = z/Az pour 
differentes demi-longueurs relatives = Az/r, la constante relative de cou- 
plage etant fixee (k = J/kc = 1/16). 

• Le cas degenere (A^ = oo) correspond au cylindre magnetique souple in- 
fini : le soliton magnetique n'etant pas confine pent s'etendre sans restriction 
aucune tout le long de I'axe de revolution et son energie magnetique ^mag 
atteindre son energie minimum minimorum topologique Eg. 

• Les autres cas (A^ < oo) illustrent le mecanisme de liberation geometrique : 
le soliton magnetique confine ne pouvant pas s'epanouir suivant I'axe de rota- 
tion refoule son exces d'energie magnetique AE^a^ = E^ag — Eq en retrecis- 
sant radialement la membrane magnetique souple ; I'energie du soliton etant 
localisee en sa vague et I'energie elastique isotrope, les deformations sont plus 
prononcees en ses branches. 
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3 Membranes de topologie torique 

3.1 Le tore : une degenerescence solitonique 

La representation la plus naturelle du tore est donnee en coordonnees 
cylindriques {p, ^, z) : 

p = i? + r cos </?, z — rsaup, (2.80) 

on le rayon de revolution R et le rayon axiale r verifient < r < R, Tan- 
gle ip variant de — tt a tt. Toutefois, il s'avere plus commode d'utiliser la 
representation suivante [Zc90] 

a sinh b a sin ri 

P = —ri . ^^—^ — — ' 2.81 

cosh — cos 7] cosh o — cos 7] 

on les nouveaux parametres constants a et b sont des reels positifs, tandis 
que V angle poloidale rj varie de — tt a tt. Les relations 

p 

{R + r){R-r) et cosh6=- (2.82) 



fournissent une interpretation geometrique immediate des nouveaux parametres 
appelons a le rayon geometrique et b Vangle d'excentricite. Reciproquement, 
les parametres naturels R et r sont donnes par 

^ 0, a ,^ 

i?= — et r = — — . 2.83 

tanh b smh b 

Enfin, la transformation entraine la relation 

tan i77 = tanh ife tan (2.84a) 



La metrique de la membrane magnetique en coordonnees peri-polaires (^, rj) 
s'ecrit alors 



g — 2 

(cosh b — cos T]) 

ainsi g^'^ — g^^ = et nous avons 



sinh^ b d^^d^ + drj^drj] , (2.85) 



g^^Vg = ^"'V^ = sinh b. (2.86) 
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Pour le tore rigide I'expression (2.4) du Hamiltonien magnetique s'ecrit done 
litteralement 



02 

— f- + sinh b sin^ 6 $^ 
sinh b ' 



+ 



. , , ^2 sin^ ^2 
smh bQl + — — $1 
' smh b ^ 



(2.87) 



Un ehangement d'eehelle approprie suivant I'axe des z donne 



= J J dc/dr/ 

-AC -T 



+ sin^ e $2 



+ 



+ sin^ e $2 



avee 



C = sinh b $, et AC = sinh b n. 



(2. 



(2.89) 



De eet Hamiltonien (2. 
rigide magnetique : 



nous deduisons les equations auto-duales du tore 



c 



e sin $ 



-s sin <l>^. 

Sans effort nous ecrivons leur versions algebriques 



(2.90a) 
(2.90b) 

(2.91a) 
(2.91b) 



En vertu du critere de Cauchy-Riemann, W est une fonction meromorphe 
[WW52] de la variable complexe cu = + e\r], le signe e etant le signe de la 
charge topologique Q. Introduisons le reseau L(2AC, \2n) = {2mA( + \2nn \ 
m,n e Z} et notons que le tore est homeomorphe a C/L(2AC, i27r) : W 
est done une fonction elliptique [WW52, Law89]. Par consequent la charge 
topologique Q de la configuration W correspond a Vordre de la fonction 
elliptique W [WW52], an signe pres. Or I'ordre d'une fonction elliptique est 
toujours superieur on egal a deux ; nous avons done 



IQI ^ 2 



(2.92) 



en accord avec le theoreme de Eells-Wood [EL78]. En outre, la representa- 
tion des fonctions elliptiques a I'aide des a-fonctions de Weierstrass [WW52, 
Law89] autorise une ecriture formelle des configurations topologiques auto- 
duales : 



W{uj) 



tan 100 e-" R^^, 

„=1<7(^-Pn) 



(2.93a) 
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ou les zeros Zn et les poles p„ doivent verifier la regie de selection 

IQI IQI 

J2^n^J2Pn- (2-93b) 

n=l n=l 

Dans la foule, ecrivons leur densite Hamiltonienne magnetique Ti-mag a I'aide 
des (-fonctions de Weierstrass : 

1 2 



H 



mag 




IQI IQI 

n=l n=l 



(2.93c) 



3.2 Sections de tore : deformation solitonique 

Restreignons nous maintenant a une membrane magnetique souple com- 
pactifiee M dont la forme spontanee est une section de tore. 

En coordonnees peri-polaires (^, rj) precedemment introduites, une telle 
membrane rigidifiee est decrite en restreignant Tangle de rotation : 

-AC ^ e ^ +Ae, (2.94a) 

ou le demi-angle doit verifier 

< A^ < TT. (2.94b) 

Encore une fois nous nous bornerons a etudier les configurations symetriques : 

= = 0; (2.95) 

le Hamiltonien magnetique (2.4) prend alors la forme attendue 

(2.96) 



avec 



-AC -TT 



C = sinh b ^ et A^ = sinh b A^. 



(2.97) 
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Enfin, la section de tore est compactifiee en imposant une densite Hamiltoni- 
enne nulle aux bords : 



Gc (C = ±AC) = 0, 

e(c = -AC) = o et e(c = +AC) 



Oh] 



(2.98a) 
(2.98b) 



II saute aux yeux que le magnetisme des sections de tore magnetiques 
souples rigidifiees est similaire a celui des sections de cylindre magnetiques 
souples rigidifiees etudie dans la Sous-Section 2.3 : seule la forme spontanee 
des membranes souples difi^ere. Autrement dit, les geometries ne sont pas sim- 
ilaires mais le choix des metriques (2.56) et (2.85) rend le magnetisme simi- 
laire. Par contre la difi^erence entre les formes spontanee est intrinseque : les 
courbures principales sont globales {i.e. constantes) pour les sections de cylin- 
dre, locales pour les sections de tore. II en decoule des parametres de forme 
pertinents de nature differente. Pour les sections de cylindre, le parametre 
de forme pertinent est le rayon r : une distance. Pour les sections de tore, 
le parametre de forme pertinent est V angle d'excentricite b : un rapport de 
distance. Nous nous attendons done a trouver des lois de deformation simi- 
laires mais aux consequences geometriques dissemblables : nous allons mon- 
trer qu'elles sont semblables. 

D'apres (2.97) et par analogic, Tangle d'excentricite b (le parametre de 
forme pertinent) doit etre libere en fonction de I'angle ^ (la coordonnee 
curviligne pertinente) : 



b{0 = bo + A{0, 



(2.99) 



avec bo V angle d'excentricite spontane et A la deformation locale. La metrique 
reste orthogonale et les relations (2.86) se lisent 



1 



y'sinh^6-F A|' 



g"^"^ = ^Jsmh^b + ^l (2.100) 



De long calculs fastidieux conduisent, en utilisant des outils classiques de 
geometric differentielle [Str61], a I'expression suivante de la courbure moyenne 
locale : 



sinh b 



^/sinhV+A| . 
+ 



sinh b + 

cosh b 



cosh b cos 1] — 1 



sinh b 



COST] 



sinh^ b + 



sinh 6A^^ -cosh 6A| . (2.101) 
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Clairement, pour des deformations suffisamment lisses, le second terme s'e- 
vanouit et le facteur d'echelle du premier terme tend vers un : I'expression 
qui s'en degage correspond alors a la courbure moyenne rigidifi,ee. En con- 
sequence de quoi, nous essayons comme courbure moyenne spontanee non 
nulle Hq la fonction 



2a Hq 7]) = sinh bo + 
+ 



cosh bo cos r] — 1 
sinh bo 
cosh b — cos 1] 



sinh^ 6 + A| 



X2 sinh bA^^ - xi cosh 6A| , (2.102) 



avec Xi et X2 des constantes phenomenologiques . Pour parachever notre mod- 
elisation, developpons la densite du Hamiltonien elastique (2.34) jusqu'a I'or- 
dre deux en A, et A^^ : 



-AC 



^ei = j dC [csA^ + 2C1X2AJ + xl^lc 

-AC 



(2.103a) 



avec 



X2 = 
Cl = 

C2 = 



1 -X2, 

1 + sinh bo cosh bo 



sinh^ bo ' 
1 + sinh bo cosh bo + sinh^ bo 
sinh^ 6n 



(2.103b) 
(2.103c) 

(2.103d) 



Precisons le changement d'echelle effectue : 



C = sinh6oC et AC = sinh6oA^- 



(2.104) 



Les deformations spontanees A qui extremisent cet Hamiltonien elastique 
(2.103) satisfont aux equations d'Euler-Lagrange suivantes [GH96] 



C2A - 2ciX2Acc + X2ACCCC = 0- 



(2.105) 



Pour interdire toutes deformations spontanees nous imposerons 



X2 = 0. 



(2.106) 
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Finalement, pour des deformations continues de faible amplitude, le Hamil- 
tonien elastique (2.34) s'ecrit 



-AC 

+AC 

= ^nkcC2 JdC{b- bof. 



(2.107a) 
(2.107b) 



Les deformations A subies par les sections de tore souples magnetiques 
et induites par les solitons magnetiques frustres geles extremisent ainsi le 
Hamiltonien total 



+AC 



= 27r JqJ / dC 



-AC 



(2 + m) - m coth 6o A + ^A^ 



- (1 + m) 2 + coth^ bo A^ sn^ (q^C | 1 + ^) 

( 1 ~ l + m C2 \ . 2 
V smh^'fto 4Kq2y 



Les equations d'Euler- Lagrange prennent alors la forme escomptee 
Agg + (1 + m)A — mB sn^ {q\ 'fn) A. — \fm j, 



(2.108) 



ou 



m — 1 + m, 
1 



mq 



1 + 



1 + m ^ C2 



2 + m \sinh bo req^^ 



B = 2 



coth^ foo 
mq2 ^ 

coth 6o 



m 



(2.109a) 

(2.109b) 
(2.109c) 

(2.109d) 
(2.109e) 



La nouvelle variable verifiant ici g = q^C- Comma pressenti les lois de defor- 
mation des sections de tore sont similaires a celles des sections de cylindre. 
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Figure 2.6: La dilatation relative At d'une section magnetique souple du 
tore de Clijford (Figure 2.3) en presence d'un 7r-soliton magnetique : la 
dilatation relative definie par (2.111) est representee en fonction de Tangle 
de rotation ^ pour differentes valeurs de la constants relative de couplage 
K = J/kc, le demi-angle de rotation etant fixee (A,^ = 7r/2). 

• Le cas degenere (k = 0) correspond a la section souple rigidifiee {kc ^ J). 

• Les autres cas (k ~ 1) illustrent le mecanisme de liberation geometrique : le 
soliton magnetique fortement confine (sinhfoo = 1) refoule son exces d'energie 
magnetique AE^ag = ^mag — Eq en retrecissant globalement la membrane 
magnetique souple ; I'energie du soliton etant localisee en sa vague et I'energie 
elastique isotrope, la contraction est moins importante en sa vague et la 
membrane souple I'epouse d'autant plus que la constante relative de couplage 
K est grande. 
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A fortiori les changements d'echelle (2.61) et (2.89) dictent la relation de 
similitude entre les parametres de forme pertinents des sections de cylindre 
(le rayon r) et des sections de tore (I'angle d'excentricite b) ; en utilisant une 
notation evidente, nous ecrivons 

En vertu de cette similitude geometrique (2.110). le mecanisme de deformation 
des sections de tore se deduit litteralement de celui des section de cylindre : la 
competition entre I'energie magnetique et I'energie elastique entraine un ac- 
croissement de Tangle d'excentricite b, ce processus s'amplifiant aux branches 
solitoniques. Observant que le produit du rayon exterieur i? = it! + r et du 
rayon interieur R = R — r reste constant {RR — selon (2.82)) definissons 
la dilatation relative du tore Xt par 

R f RV^ 

oil Ro et Rq sont associees a la forme spontanee. Les relations de transforma- 
tion (2.83) permettent d'etablir une relation simple entre la dilatation relative 
du tore Xt et V angle d'excentricite b : 

tanh ifon 

At = r^. 2.112 

tanhift ^ ^ 

3.3 Deformations symetriques du tore de Clifford 

Le tore de Clifford magnetique souple serait assouvi puisqu'il saturerait 
son energie elastique selon la conjecture de Clifford (Figure 2.3) tout en 
saturant son energie magnetique du fait de sa degenerescence solitonique 
(Sous-Section 3.1). Par deformations symetriques du tore de Clifford nous 
entendrons les deformations subies par le tore de Clifford magnetique souple 
et induites par un soliton magnetique symetrique : en fait il s'agit du cas 
degenere = tt d'une section du tore de Clifford vu dans la Sous-Section 
precedente. La Figure 2.7 illustre cette approche. 
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Figure 2.7: Deformation symetrique du tore de Clifford magnetique souple 
en presence d'un 27r-soliton magnetique symetrique : I'angle d'excentricite bo 
du tore spontane (representee par la grille exterieure) verifie sinh bo = 1 par 
definition (FIGURE 2.3), la constante relative de couplage k= J /kc est egale 
a 2, enfin la deformation locale A est magnifiee par un facteur 10. 
Le soliton frustre se libere de sa symetrie imposee en contractant globale- 
ment la membrane magnetique ; les deux renflements diametralement op- 
poses correspondent aux deux vagues du soliton : son energie magnetique 
Emag s'y conccntraut independamment de la geometric de la membrane, il 
refoule I'energie magnetique accumulee en ses branches. 
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La premiere partie esquisse la construction de nouveaux degres de liberte an- 
tiferronnagnetiques pour un reseau regulier de basse dinnension. Des polytopes 
reguliers du reseau ils heritent de la hierarchie des algebres a divisions : le nou- 
veau parametre d'ordre local suit la hierarchie algebrique satisfaite par le modele 
cr non-lineaire. Le reseau de I'etat antiferromagnetique (le « sur-reseau ») appa- 
raTt alors comme un reseau de noeuds (les « sur-nceuds ») avec liens (le champs 
potentiel vecteur associe au champs de jauge locale). Une etude approfondie 
du nouveau parametre d'ordre local pris a la limite continue semble un prelim- 
inaire necessaire au passage a une theorie de champs continue. Notons enfin 
que I'argument naTf qui consiste en dimension deux (algebre des complexes C) 
a envoyer le plan-temps compactifie 5*^ sur la sphere unite pour degager 
invariant topologique de Hopf n est plus pertinent en dimension trois (algebre 
des quaternions H) puisque la hierarchie algebrique impose alors d'envoyer la 
sphere unite 5"^ sur la sphere unite S'^ [EHH+90, Nak90] : une etude du reseau 
cubique antiferromagnetique devrait permettre d'afFiner la prescription esquis- 
see, d'autant plus que reseau cubique antiferromagnetique peut-etre soit assouvi 
soit frustre [Kaw98]. En dimension quatre (algebre des octonions O), I'invariant 
topologique de Hopf decompte les recouvrements de la sphere unite -S"^ lorsque la 
sphere unite S^^ est parcourue : une etude du reseau 4-cubique et 4-triangulaire 
[Cox73] devrait permettre de celer la prescription et une meilleur comprehension 
du mecanisme algebrique sous- jasent puisque I'algebre des octonions O n'est 
ni associative ni commutative mais seulement alternative — I'alternativite etant 
commune aux quatre algebres reelles a divisions A<i. 

La seconde partie s'evertue a decrire la membrane magnetique comme un 
systeme physique possedant une double obstruction topologique. La frustra- 
tion geometrique s'interprete alors comme la competition entre les deux ordres 
topologiques. En outre cette approche suggere la construction d'un Hamiltonien 
elastique qui admettrait une decomposition «a la Bogomol'nyi*, c'est-a-dire 
qui leverait la classification topologique des surfaces comme le modele a non- 
lineaire (le Hamiltonien magnetique) leve la classification des configurations mag- 
netiques. Le Hamiltonien heuristique de Helfrich-Willmore qui mimique une telle 
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decomposition a permis de mieux saisir le mecanisme de liberation geometrique. 
Les deformations d'une membrane souple ont ete apprehendees en gelant une 
configuration magnetique de ladite membrane souple rigidifiee : un retrecisse- 
ment global avec des renflements localises aux retournements des spins a ete mis 
en evidence pour les membranes magnetiques frustrees. Ce resultat demande a 
etre confirme en elaborant un Hamiltonien elastique independant de la topolo- 
gie de la membrane souple (c'est-a-dire «a la Bogomol'nyi ») et en abordant la 
competition entre les deux ordres topologiques conjointement. Notons que I'ex- 
istence d'un tel Hamiltonien elastique souleverait I'existence d'Hamiltoniens «a 
la Bogomol'nyi ». 

Le processus microscopique qui se degage de la prescription suggeree dans 
la premiere partie et la competition topologique introduite dans la deuxieme 
partie elargissent le champs d'investigation initial en fournissant une origine mi- 
croscopique a la transmutation de Fermi-Bose et un moyen de traquer la physique 
exotique etudiee lors de cette these. Des simulations numeriques du modele de 
Heisenberg sur des reseaux bidimensionnels souples spontanement courbes de- 
vrait ainsi exhiber des distributions topologiques de spins dont I'ordonnancement 
local permettrait de preciser le processus microscopique. De meme I'etude des 
transitions de phase entre les difFerents etats antiferromagnetiques topologiques 
devrait egalement s'averer utile. 



Annexe A 

L'equation de sinus-Gordon 



Cet Annexe presente des solutions periodiques de l'equation double de sinus- 
Gordon(DSG) obtenues elegamment a partir de solutions periodiques de l'equa- 
tion simple de sinus-Gordon (sg). 



1 L'equation simple de sinus- Gordon 

Considerons Vequation de sinus-Gordon (sg) 

©p^ = sin©cos©. (A.l) 

En integrant une premiere fois l'equation differentielle du second degres 
(A.l), nous lisons 

©J = sin^©-Fm me[0,+oo[, (A. 2a) 

ou le parametre m est une constante d'integration. II est alors aise d'ecrire 
l'equation satisfaite par la fonction reciproque Q — Q (©) : 

^ - , " . ^ - ±1- (A-2b) 

d© Vm + sin" ^ ' 

Notons a (■ I m) la fonction croissante solution de (A.l) caracterisee par le 
'parametre m et fixee a I'origine, c'est-a-dire verifiant a (0 | m) = |. La fonc- 
tion reciproque de a (■ | m) est donnee par Vintegrale elliptique 



a 



-'{e\m)^ f^=±==, (A.2c) 
' 'm -|- sin (f 



2 

solution immediate de (A. 2b). Ainsi les solutions SG envisagees s'ecriront- 
elles, avec des notations naturelles, 

e (g) ^ea{Q\m)+ c. (A. 2d) 
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En observant que 

(6' + TT I m) + 2a~^ (0 | m) = {0 \m) , 
nous definissons la quasi quart-periode par Vintegrale elliptique complete 

TT 

K^{m)= l^=±==, (A.3) 

m + sin (/? 



qui decroit strictement de oo vers en fonction du parametre m. La fonction 
a{- \m) verifie done la relation 

a{g + 2qKa \ m) — a {g \ m) + qn q € (A. 4a) 

de plus, 

a{-Ka)=Q, «(0) = f, a{Ka) =7r; (A.4b) 

enfin, par construction, 

I — q; {—g \ m) — a {g \ m) — ^. (A. 4c) 

Clairement le parametre m controle la quasi periodicite des solutions. Lorsque 
le parametre m tend vers 0, cette quasi periodicite tend vers I'infini, et la 
solution a{- \rn) devient un soliton localise a I'origine. 

Uequation de sinus-Gordon a ete resolue sans effort aucun, mais les re- 
lations ainsi obtenues sont peu commodes a utiliser en calcul symbolique ou 
numerique. Des expressions explicites peuvent etre ecrites en termes dlntegrales 
canon iques elliptiques et des fonctions elliptiques de Jacobi [WW52, Law89]. 

En effet, en effectuant le cliangement de variable judicieux [Sut95] 
sin© = dn (ti I 1 + m), l'equation difTerentielle originelle (A.l) prend la forme 
simple v?Q — 1. Immediatement, nous lisons 

sina I m) = dn I 1 + m) . (A. 5) 

En substituant cette relation dans (A. 2a), nous obtenons 

o;^ {q \m) = {1 + m) ct? {q\ 1 + m) . (A. 6) 

La trigonometric classique et jacobienne permettent de calculer la variante 
cosinus de (A. 5) : 



cos a {g\ m) — — vT+m sn(f) | 1 + m) , (A. 7a) 

= -sn (^Vl + m I 1/1 + m) . (A. 7b) 



2 L 'equation double de sinus- Gordon 
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Pour finir, ecrivons la variante en tangente de I'angle moitie : 

tan ia I m) = nd I 1 + m) + Vl + rra sd | 1 + m) , (A. 8) 

qui prend la forme degeneree, quand m tend vers 0, tan {g \ 0) ~ e^. 

Concernant la quasi quart-periode (m) , le rapprochement des relations 
(A. 4a) et (A. 5) permet de I'identifier avec la quart-periode K (1 + m) des 
fonctions elliptiques de Jacobi pq(- | 1 + m). Precisement, nous avons 

Ka (m) = K (1 + m) , (A.9a) 

= -r^=^ K (-^) . (A.9b) 
VTTm \l + mj ^ ^ 

En calcul symbolique, K correspond a Vintegrale elliptique de premiere espece. 



2 L 'equation double de sinus- Gordon 



Uequation double de sinus-Gordon (dsg) s'obtient en ajoutant un terme 
de couplage a Y equation simple de sinus-Gordon : 

= sine cos© + 7sine. (A. 10) 

Le changement de variable 

^ = tani© ^e[0,+oo[ (A.ll) 

transforme I'equation differentielle non lineaire (A. 10) sous la forme plus 
commode 



( 



1 + r t 



:2t 



^« + l(l + 7)-|(l-7) t' 



(A.12) 



Cette equation differentielle admet deux solutions triviales de type soliton 



le double soliton t — cosh -^g 



I) et le soliton t 



(7 = 0). 



Clairement, la transformation {g, t) — > (Xg^iit) et 1' absence de couplage 
7 ne modifient pas la nature de la nouvelle equation DSG (A.12). Cette re- 
marque suggere fortement d'injecter dans (A.12) la fonction essai 

t — tan icKo tan ia {g/vQ \ m). (A. 13) 

De cette manipulation fastidieuse se degagent un systeme d'equations dont 
s'extraient les expressions de Vq et de tan \aQ : 



tan itto- 



\ 



2 + m - ^Jm'^+^ (1 + m) 7^ 



2(l-f 



14-^ 27 -|- m — m^-|-4 (1 + m) 7^ 
1 - 7 27 - m + ^m2+4(l + m) 7^' 



(A.14a) 
(A.14b) 
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Lorsque le couplage 7 disparait, la longueur vq et tan iao tendent vers I'unite : 
sans surprise la fonction essai (A. 13) approche la fonction SG associee. Ainsi 
ces solutions periodiques se revelent-elles etre une generalisation des solutions 
a (• I m). Notons a (• | m, 7) les solutions (A. 13) ; la formule de Tangle moitie 
correspondante s'ecrit 

tan ia (g | m, 7) = tan icto nd (g/vo \ 1 + m) + Vl + rn sd (g/vo | 1 + m) , 

(A.15) 

ou Tangle a I'origine ao et la longueur solitonique vq verifient (A. 14). En 
consequence de quoi, la quasi quart-periode generalisee (rn, 7) satisfait 

Ka {m, 7) = vo Ka (m) . (A. 16) 

La relation (A. 4a) se generalise alors comme suit : 

a{g + 2qKa\m,^)=a{g\m,{-iy^) + qTT q e Z; (A. 17a) 

de plus, les identites (A. 4b) deviennent 

a{-Ka)^0, a{0)^ao, a{Ka)^n; (A. 17b) 

finalement, la relation de parite (A. 4c) prend la forme 

I — q; {—g \rn,^) — a{g\m, —7) — |. (A. 17c) 



3 Evaluation numerique du parametre 

En physique les conditions aux limites sont generalement fixees. Ce qui 
revient ici a se donner la quasi quart-periode K^. Or il n'existe pas dans 
la litterature de relation ou d'algorithme numerique permettant d'evaluer la 
fonction reciproque : une methode numerique ad hoc doit etre suggeree. 

Pour commencer, notons que Vintegrale elliptique de premiere espece K (m) 
pent s'exprimer sous la forme d'une fonction hyper geometrique [AS72, WW52] : 

KW = fF(ii;i;m). 



3 Evaluation numerique du parametre 
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Cette remarque et la relation (A. 9b) incitent a ecrire comme une serie 
entiere du module reduit x = l/Vl + ttT'- Ces considerations conduisent a 
I'expression 



ou les symboles de Pochhammer [AS72] sont utilises. Cette serie entiere est 
alors formellement inversee pour obtenir la serie entiere representant la fonc- 
tion X (Ka). Un scripte (code en Maple V) permettant de calculer les premiers 
coefficients de cette serie s'ecrirait comme suit : 

restart ; 
Ninax:=100; 

a:=(n)->piecewise(type(n,odd) , 

(pochhammer ( 1/2 , iquo (n , 2) ) /pochhammer ( 1 , iquo (n , 2) ) ) "2 , 
type (n, even) ,0, 
'a'(n)) : 

powseries [powcreate] (K(n)=a(n) ,K(0)=0) : 
M:=powseries [reversion] (K) : 

SineGordonM [ApproximationOrder] : =Nmax : 
SineGordonM [nestedPowerSerie] :=table() : 
for n from to Nmax+1 do 

SineGordonM [nestedPowerSerie] [n] :=M(2*n+l) : 
od: 

saveCSineGordonM, "SineGordonM. m") : 

De cette fagon, la relation (A.18) permet V evaluation numerique de la fonc- 
tion reciproque K^ ^ pour des valeurs numeriques suffisamment petites. 
Pour des valeurs trop grandes, les inegalites suivantes [AVV90] 




(A.18) 



In 



\/l — m 



^ K(m) ^ In 




[0<m<l] 



(A.19) 



s'averent utiles (et necessaires) pour fournir a notre procedure f solve fa- 
vorite un intervalle initial suffisamment fin pour resoudre (avec succes) I'equa- 
tion Ka = X K (x^) numeriquement. 




Figure A.l: La fonction reciproque de la quasi quart-periode de sinus- 
Gordon Kq.. 



Annexe B 

L 'equation inhomogene de Lame 



Cet Annexe se contente de presenter une solution particuliere de I'equation 
inhonnogene de Lanne(lL) avec un terme inhonnogene constant. La construction 
de cette solution particuliere s'inspire de la construction des fonctions dites de 
Lame [Erd55, Ars64]. 



1 L 'equation aux differences 

Soit Vequation inhomogene de Lame (il) 



ou les amplitudes A et B sent des parametres constants. A priori, le terme 
inhomogene j est un reel non-nul et le parametre m un reel positif. 

Comme la transformation reelle de Jacobi (g, m) — > {Q^/m, m~^) laisse 
I'equation differentielle lineaire du second ordre (B.l) inchangee, nous pou- 
vons suppose m compris entre et 1. En outre, puisque la fonction sn^ {g \ m) 
est paire, I'equation differentielle (B.l) reste inchangee sous la transformation 
de parite g—>-—g, done une solution particuliere pertinente sera une fonction 
paire en g. De maniere analogue, la parite en ^ — K (m) de sn^ {g \ m) invite 
a chercher une solution particuliere egalement paire en ^ — K(m). Ces re- 
marques preliminaires faites, et tenant compte de la nature meme de (B.l), 
cherchons une solution particuliere sous la forme 



Kqq + (1 + m)A — mB sn^ {g \ m) A — ^/m j, 



(B.l) 



oo 



^ On sn^" {g\ m) . 



(B.2) 



n=0 
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La substitution de la fonction formelle (B.2) dans (B.l) impose a la suite 
des coefficients reels a„ la relation de recurrence a trois termes [Erd55, Ars64, 
Gau67, LW92, PTVF92] suivante : 

- A 

Oa+l - (1 + m) an 
[2n + 2)[2n + 1) 

(2n - 2)(2n - 1) - B ^ , 

+ " (2« + 2)(2n + l) ^-' = ° <^-'''' 
pour n = 1,2, . . . ; la relation d'amorgage etant 

2ai + (1 + m)A oq = a/^ J- (B.3b) 

L'equation caracteristique de (B.3) se lit immediatement : 

- (1 + m) a + m = 0; (B.4) 

elle admet deux racines evidentes, m et 1. Done d'apres le theoreme de Perron 
[Gau67, PTVF92], pour ^ m < 1 (i.e. m ^ 1), l'equation aux differences 
lineaire du second ordre (B.3) admet deux solutions lineairement indepen- 
dantes, a savoir la solution minimale et la solution dominante, dont les rapports 
On+i/ On tendent, respectivement, vers m et 1. 



2 La solution minimale 

Ecrivons la relation de recurrence (B.3a) sous la forme 

(2n-2)(2n-l)-B 
_ ""^ (2n+2)(2n+l) .p. 

On-l -(1+m) (2„+2)(2n+l) + ^+1/ ^ 

Si la nouvelle relation de recurrence (B.5) est reitere indefiniment a partir 
d'un certain rang n ^ 1, une suite de fractions est alors engendree [LW92]. 
Cette suite, selon le theoreme de Pincherle [Gau67, LW92, PTVF92], con- 
verge vers le rapport On/ ftn-i ou les coefficients et On-i sont effectivement 
ceux de la solution minimale. Notons r„ la fraction continue precedente, les 
coefficients On sont alors donnes par 



an=ao]Jrp [n=l,2, ...], (B.6a) 

p=i 



avec, selon (B.3b), 



«o = 7;^ tt; (B.6b) 

2ri + {1 + m)A' ^ ' 
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de plus, avec des notations usuelles [LW92], le rapport r„ s'ecrit 



K 



(2p-2)(2p-l)-B 
(2p+2)(2p+l) 



p=n \ — (1+to) 



[n = l,2,...] 



(B.6c) 



(2p+2)(2p+l) , 



Nous adopterons la notation L (■ | m\A,B,i) pour designer la solution 
particuliere de (B.l) associee a la solution minimale de I'equation aux dif- 
ferences (B.3). La relation (B.6a) permet d'ecrire cette solution particuliere 
sous la forme de Horner 



L {q I m;A,B,j) = oq 



1 + ri sn^ (g \ m) 



1 + r2 sn^ [q I m) [• 



qui doit etre penser comme une « forme de Horner continue » . Nous ecrirons 
done, avec des notations evidentes. 



oo 

L I m; A, B, j) = Oo H [l + sn^ (q I rn) 



n=l 



(B.7) 



L'expression (B.7) n'est vraie que pour < m < 1. En fait, I'invariance de 
y equation inhomogene de Lame (B.l) par la transformation reelle de Jacobi : 



L(q\ m;A,B,j) ^L(^Q^/m\ m 



(B.8) 



3 Evaluation numerique 

La construction de L (^ | m]A,B,j) se prete naturellement a I'elabora- 
tion d'un algorithme numerique. Cependant, revaluation de la fraction con- 
tinue (B.6c), autrement dit des rapports r„ dans (B.7), apparait clairement 
comme la principale difficulte : V algorithme de Lentz dans sa version modifiee 
[PTVF92, Len76, TB86, CJT93] doit etre suggere. 
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